Impurity in a d-wave superconductor: Kondo effect and STM spectra 
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We present a theory for recent STM studies of Zn impurities in the superconductor 
Bi2Sr2CaCu208+5, using insights from NMR experiments which show that there is a net S = 1/2 
moment on the Cu ions near the Zn. We argue that the Kondo spin dynamics of this moment is 
the origin of the low bias peak in the differential conductance, rather than a resonance in a purely 
potential scattering model. The spatial and energy dependence of the STM spectra of our model 
can also fit the experiments. 
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Recent progress [nlM] in scanning tunneling microscopy 
(STM) of surfaces of the high temperature supercon- 
ductor Bi2Sr2CaCu208+5 has given us a high resolution 
probe of its electronic correlations. Especially notable 
have been studies [y in the vicinity of Zn impurities on 
the Cu sites, and the site-specific information on the vari- 
ation in the pairing correlations. The Zn^^ ion has spin 
5 = 0, and so it appears natural to interpret ||l|,|[ the ex- 
periments using a theoretical model ^;Q| in which each 
Zn site acts a potential scatterer of Cu 3d electrons form- 
ing a d-wave BCS state. However, several key experimen- 
tal features do not fit easily in this simple picture. There 
is a large peak in the differential conductance close to 
zero bias near the Zn site: it is tempting to identify this 
peak with a quasi-bound state in the potential scatter- 
ing model ||], but such a state appears at low energies 
only for a range of very large potential values depending 
upon microscopic details &fl]. Furthermore, the spatial 
dependence of the zero bias peak is unexpected |g]: it 
is largest on the Zn site, and with a local maximum on 
the second neighbor Cu sites (see Fig 3c in Ref |^), and 
this differs from the general expectations of the potential 
scattering model [0,0 . Finally, the observed spatially in- 
tegrated spectrum is asymmetric between positive and 
negative bias ||l|], while the potential model predicts ap- 
proximate symmetry |9[| . 

We shall address these questions here using a rather 
different theoretical model which was originally moti- 
vated by other experiments on Zn or Li impurities (the 
Li"*" ion is also S = 0). A series of beautiful NMR ex- 
periments |10|-[l5|] have clearly shown that each impurity, 
despite having no on-site spin, induces a local, unpaired 
S = 1/2 moment on the Cu ions in its vicinity at interme- 
diate energy scales. In the underdoped regime, this can 
be understood [ [lO|Jl^ by the confining property of the 
host antiferromagnet, in which the impurity is a localized 
"holon" which binds the moment of a 5 = 1/2 "spinon" . 
At larger doping, a related picture can be developed by 
analogy with the theory p7| for moment formation in 
the disordered metallic state of Si:P — small variations in 
the potential combine with strong local correlations to 



induce very localized spin excitations. Other theoret- 
ical perspectives [|l8| , |l9[ on local moment formation in 
the cuprates have also been given. Further evidence for 
S — 1/2 local moments near Zn sites appears from neu- 
tron scattering experiments |Q : it has been argued |2^] 
that these are required to explain the strong effects of 
a small concentration of Zn impurities on a "resonance 
peak" in the spin dynamic structure factor. 

This paper will study a 5 = 1/2 local moment near 
the Zn/Li site coupled by exchange interactions to the 
fermionic 5 = 1/2 excitations of a d-wave superconduc- 
tor. We shall describe the effects of the Kondo screening 
of the local moment, by the fermionic excitations, on the 
STM spectra: we shall show that our model leads natu- 
rally to a peak in the differential conductance at a low 
bias of order the Kondo temperature, with a spatial dis- 
tribution which can fit the experiments. Our spatially 
integrated spectrum also has a bias asymmetry which in- 
creases with decreasing doping. 

Wc begin by describing our model Hamiltonian, H ~ 
Hbcs + -ffimp, for a single Zn impurity ||2l|. The first 
term describes the host superconductor, which we model 
by a simple BCS Hamiltonian 

i^Bcs - ^ *I [(^k ~ m)t" + AkT^] *k. 

k 

Here ^k = (cktiC_k|) i^ ^ Nambu spinor at momentum 
k = {kx, ky) (cka annihilates an electron with spin a on 
a 3d orbital), r^'^'^ are Pauli matrices in particle-hole 
space, and /i is a chemical potential. For the kinetic 
energy, Ek we have first (t), second (t'), and third {t") 
neighbor hopping, while we assume a d-wave form for 
the BCS pairing function Ak = (Ao/2)(cosfc2: — cosfcj^). 
The Zn impurity is at r = Tq, and is described by 



i?imp = E Kir)S-ciir)^c^ir) + C/4(ro)c„(ro), 

rGA/' 



where ^/N^Ca{r) = X^fe '^kaC*'''" {Ng is the number of 
sites in the lattice) annihilates an electron at site r, a are 
the Pauli matrices in spin space, and A/" is a set of sites 



in the neighborhood of rg. The spin degree of freedom 
induced by the impurity is represented by the S = 1/2 
operator S. The on-site potential scattering of the Zn ion 
is represented by U: previous analyses of STM spectra 
of Zn ions |^,^^ included only this second term in Hi^p 
and omitted [9[ the degrees of freedom represented by 
S. We assume that effects due to spatial variations in 
the self-consistent pairing amplitude near the impurity, 
along with those due to Hartree renormalizations from 
the Coulomb interactions, have been absorbed into the 
effective parameters K{r) and U (as in ||l^); therefore, as 
the dopant charge carriers are moving in a background of 
ions that are both doubly charged (Cu"'"*' and Zn^^), we 
expect U to be negative and measuring mainly the shift 
in the d level energy on the Zn (along the same lines, U 
should be larger and positive for Li^ impurities). 

An important ingredient in our computation, which in- 
fluences the spatial form of the STM spectrum, is the r 
dependence of the spin-dependent interaction K{r). This 
is quite difficult to determine from ffist principles, and 
we will instead use the spatial information obtained by 
analyses of Knight shifts in NMR experiments ||l4| , [l5| . 
The S = 1/2 moment is found to be concentrated on the 
4 Cu nearest neighbors of the Zn ion (see Fig Ig in fl^ ] 
and Fig 19 in ||lj]), and the moment on the Zn site itself, 
r = rg, is negligible. We expect that the dopant holes will 
have a strong preference to reside on a Cu site with a mo- 
ment (as the Cu spins on the other sites have paired with 
each other and gained exchange energy); this attraction 
is realized by K{r), and we take K{r) = Ki for the 4 sites 
r — ro = (±1,0), (0,±1), while K{ro) — 0. In principle, 
it is not difficult to also include the smaller K{r) values 
at larger values of |r — ro|, but we will neglect them here 
for simplicity. We also note that while the above form 
for K{r) is reasonable, there is no justification for a cor- 
responding form for the potential scattering term, which 
is surely largest at r = ro, as in Tfimp- 

The Kondo effect in H is related to that in a class of 
models which have been much studied recently [p2|-p4[: 
these models have a single spin coupled to a fermionic 
bath whose local density of states, p^(e), vanishes as a 
power law near the Fermi energy, pi{e) ~ |e|'', where e is 
measured from the Fermi level. The usual Kondo effect 
in a Fermi liquid corresponds to r = 0, while the present 
d-wave superconductor has r = 1. The value of r and 
the presence or absence of particle-hole symmetry |2J] 
play a key role in determining the low energy physics, 
and a comprehensive phase diagram has been presented 
by Ingersent and collaborators |gj|. For r — 1 and with 
perfect particle-hole symmetry (for H this corresponds 
to i' = 0, /i = 0, and C/ = 0) there is in fact no Kondo 
effect: the spin is free at low energies for all values of 
the exchange, K. However, there is a quantum phase 
transition at a finite magnitude of particle-hole symme- 
try breaking and at a finite K, to a phase in which the 



spin is Kondo screened below an energy scale Tk- The 
universal critical theory for this quantum critical point is 
not known and we shall not discuss it here. 

Here, we shall describe the dynamics of H using a 
"large iV" approach |2|,||,|2|. While this method has 
numerous artifacts near the quantum-critical point just 
mentioned, and is not quantitatively accurate, it does 
capture the qualitative physics of the Kondo screened 
phase in an effective manner, and this is our primary in- 
terest. We represent the spin S* by a fermion, fa, which 
obeys the constraint flfa = 1. We impose this by a La- 
grange multiplier A, and decouple the exchange interac- 
tions by complex Hubbard-Stratonovich fields ip^] A and 
(^r are approximated by their static, real saddle point val- 
ues. The physical quantities are expressed in terms of the 
Green's function of the Nambu spinor F = (/|, /J), which 
in Matsubara frequency, a;„, is T{lj„) — (i^(a;„)F'l'(a;„)). 
We have 
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i^r'/5r'T^G(r, r',u;„)T'= 



where G is the ^ Green's function with potential scat- 
tering alone: 

G(r,r',c^„) = G°(r - r',c^„) ~ C/G'^(r - ro,c^„) 

xr^[l + C/G"((0,0),c^„)Tl"'G°(ro-r',L^„); (1) 

G° is the Green's function of the host 7VsG°(r,a'„) — 
J2i,e'^-''[-iu;n + {sk - t^y'' + Afer^]"!. The values of A 
and (^r were obtained by numerically solving the con- 
straint equations T^^ Tr[T^T(w„)] — and ipr = 

-(Tif (r)/2) E.^ Er'eA^>r'Tr[T^T(a;„)r-G(r',r,^„)] 
(T is the temperature, and Boltzmann's constant fc^ = 
1). Finally, the tunneling density of states (DOS) was 
obtained as Im Tr(G(l + t^))/2, where G is the fuU * 
Green's function 

G(r,r',Lj„) = G(r, r',t^„) 

-f Y^ ipsVs'G{v,S,UJn)T'T{uJn)T'G{s' ,r' ,UJn). (2) 

s,s'ej\f 

We now describe our numerical results. For small Ki, 
the constraint equations have only the solution if^ = Q. 
In the present large N approach, S is completely de- 
coupled from the fermions at such a saddle point. This 
solution corresponds to the free spin phase, and the 
STM spectra of the large N theory are identical to 
that of the purely potential scattering model. To in- 
clude the spin dynamics, we take Ki > Kic below so 
that ifr 7^ 0; we found that the lowest free energy sad- 
dle points had a d-wave pattern with tp^ = +[—]ip for 
r — ro = (±1, 0)[(0, ±1)]. Provided particle-hole symme- 
try is absent, the onset of a non-zero tp corresponds to 
a phase transition to the Kondo screened phase at Kic 
(for the particle-hole symmetric case, there is no Kondo 
screening p3] , and the large TV equations show that this 



is so, even for Ki > Kic). The assumption of eventual 
low energy Kondo screening is also in accord with indi- 
cations in NMR jlj] and thermodynamic measurements 
p7[ . The large N value for Kic is believed to be an over- 
estimate: this are other shortcomings will be addressed 
in a forthcoming work employing an alternative "non- 
crossing approximation" p6| |. 

For Ki > Kic, the energy dependence of the tunnel- 
ing spectrum is dominated by the form of the scattering 
matrix T(lo); the imaginary part of T(u;) shows a pro- 
nounced maximum at an energy wq (which becomes a 
very sharp peak near the Kondo transition), so the lo- 
cation of possible peaks in the local DOS is given by 
ibwo- On the other hand, \uJo\ can be also be identi- 
fied with the Kondo temperature, Tk- To see this, con- 
sider the local impurity susceptiblity, which is given by 
Xioc - -(r/4)X;„„Tr[TK)2]. Assuming the simple 
pole structure T(ujn) — {—iujn + ojqt^)^^ it is evident 
that xioc changes its character at a T of order |wo| ~ Tk- 

We show in Fig |l| the energy and spatial dependence 
of the local DOS of H for a typical set of parameters 
at 14% hole doping. Most significant is the pronounced 
peak at rg (on the Zn site) at a bias, oj — —loq. This was 
a robust feature of our results: the peak retained a small 
negative bias and width for a wide range of Ki > Kic and 
U < 0. Away from the central site, the peaks were much 
weaker, and varied slightly depending upon U and dop- 
ing: for —0.4 eV < U < at optimal doping (and more 
robustly at lower doping) we obtained the alternating in- 
tensity pattern in Fig |l|b |28) . The spatial integral of the 
spectrum at negative bias is larger than that at positive 
bias, and the ratio is quoted in Fig n^. For comparison we 
show in Fig the analogous results for a purely poten- 
tial scattering model Q (Ki — and S absent), where 
we see dramatically different features. Now the largest 
peak is on the first neighbors 0]; notice also the change 
in the color scale in Fig 0, and that this peak is not as 
pronounced as the Zn site peak in Fig m. Further, we 
had to choose a very large value of |C/| to make the peak 
sharp and at a small negative bias. Both these features 
are sensitive to variations in the values of the doping [|6| , 
the band structure {f /t, t" /t) and U: it is not difficult 
to find broader peaks at higher energies, and to switch 
the largest peak to positive bias. 

We make a few remarks on the spatial dependencies in 
Fig y. The dominant contribution at peak energy comes 
from the second term on the right hand side of (0), and 
both the normal (t^) and anomalous (t^) components 
of G are largest for sites on opposite sublattices. Using 
the location of the four sites in Af and the d-wave struc- 
ture of both the superconducting gap function Ak and 
the ifr fields, we see that at Tq the normal (anomalous) 
Green's functions destructively (constructively) interfere 
leading to a large peak at —uq. This interference is also 
responsible for the asymmetry of the spatially integrated 
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FIG. 1. Calculated tunneling density of states for the 
Kondo impurity model H = Hbcs + J^imp at 14% hole dop- 
ing with a realistic band structure (t = 0.1 eV, i' — — f/4, 
t" = t/8), Ao = 0.04 eV, Ki = 0.21 eV, n = -0.09 eV and 
U — —0.05 eV. For these parameters Kic = 0.1 eV. (a) Tun- 
neling spectrum versus sample bias for the impurity site (red) 
and the nearest (blue) and second (green) neighbor sites, (b) 
Spatial dependence of the differential conductance at the bias 
uj = —1.9 meV {— —uJo)- The ratio of the spatial integrals at 
±u! is 1.5, and this value increases with decreasing doping. 



spectrum. (We note that for an s-wave pattern in the 
ipr fields the interference mechanism is similar, with the 
difference that the peak at Tq appears at +uJo.) An anal- 
ogous discussion [|| can be applied to Fig g, where inter- 
ference of scattering from different r is absent. We also 
note that the coherence peaks, appearing in the DOS at 
the gap energy, are much suppressed at the impurity site 
and its neighbors, although we did not include a spatial 
variation of the pairing amplitude. Furthermore, upon 
moving to the metallic state by setting Aq = 0, our ap- 
proach is closely related to that of ||2^] , and then the large 
DOS around the Fermi level leads to a Fano lineshape. 

To summarize, this paper has introduced a new model 
{H = i^BCS + ^imp) for STM spectra of Zn/Li impurities 
in the high temperature superconductors. The model ap- 
pears to be required to explain NMR [[lO|-[l^] and neutron 
scattering measurements |^,^ on the same system. It 
is therefore satisfying that our model also leads to a low 
bias peak in the STM; the measured bias of this peak 
(1.5 meV) suggests that the Kondo screening of the mo- 
ment (as can be measured directly in NMR) occurs at a 
temperature of order 15 K. 
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FIG. 2. Parameters as in Fig hi, but for a purely potential 
scattering model (//imp has A'l = and no S) with U — —4 
eV. The very large \U\ is necessary for a low bias peak. The 
ratio of the spatial integrals at ztwo is now 1.05, and this value 
remains around unity with decreasing doping. 

We conclude by noting related issues: (i) Our model 
admits a simple generalization to Ni impurities. The 
Ni++ ion has 5 = 1, and so we add an additional S = 1 
degree of freedom ^ni ■ A simple choice for a Hamiltonian 
is if -I- K'S^i ■ S. Depending upon the values of K' and 
K{r), a rich variety of behaviors appear possible. A likely 
possibility, suggested by the NMR experiments ||l^, is 
that S and S'ni combine to form a 5 = 1/2 moment; the 
Kondo coupling of this effective moment is expected to 
be ferromagnetic [Q — so a, S — 1/2 moment remains 
unscreened. The STM spectra will interpolate between 
the results described here, and those in the presence of a 
static local magnetic field Ig^l i and this appears to be the 
case [^. (iz) In a recent study [|l| of the broadening of 
the collective spin resonance mode [ pO| by Zn impurities, 
the Kondo screening of S by the fermionic quasiparticles 
was neglected. This is valid as long as the energy of the 
resonance mode. Ares ~ 40 meV, is larger than Tr-; this 
condition is obeyed by the Tk values discussed above. 

We thank H. AUoul, W. Atkinson, A. Balatsky, 
S. Davis, G. KhaliuUin, A. MacDonald, and S. Pan for 
valuable discussions, and US NSF Grant No DMR 96- 
23181 and the DFG (VO 794/1-1) for support. 

[1] E. W. Hudson et ai, Science 285, 88 (1999); S. H. Pan 

et ai, Nature 403, 746 (2000). 
[2] A. Yazdani et ai, Phys. Rev. Lett. 83, 176 (1999). 



[3] 



[4] 
[5] 
[6] 
[7] 



[8] 
[9] 



[10] 
[11] 

[12] 

[13] 
[14] 
[15] 
[16] 

[17] 

[18] 

[19] 
[20] 
[21] 

[22] 



A. V. Balatsky, M. I. Salkola, and A. Rosengren, Phys. 
Rev. B 51, 15547 (1995); M. I. Salkola, A. V. Balatsky, 
and D. J. Scalapino, Phys. Rev. Lett. 77, 1841 (1996). 
H. Tsuchiura et ai, J. Phys. Soc. Jpn. 68 2510 (1999). 
W. A. Atkinson et ai, Phys. Rev. Lett. 85, 3926 (2000). 
L Martin and A. V. Balatsky, |cond-mat/0003142. 



[23] 
[24] 



[25] 
[26] 

[27] 
[28] 



[29] 
[30] 



[31] 



If the host band structure is particle-hole symmetric (not 
the case at optimal doping), the potential scattering {U 
in -ffimp) has to be near infinity to obtain a low energy 
quasi-bound state, but in general, some finite and non- 
universal value of U is required. 
We thank H. AUoul for emphasizing this to us. 
M. E. Flatte and J. M. Byers, Phys. Rev. Lett. 78, 3761 
(1997); M. E. Flatte, Phys. Rev. B 61, 14920 (2000): this 
paper included an effective magnetic field at ro but did 
not account for the Kondo effect. 
A. M. Finkelstein et al., Physica C 168, 370 (1990). 
H. AUoul et al, Phys. Rev. Lett. 67, 3140 (1991); 
A. V. Mahajan et ai, ibid 72, 3100 (1994). 
J. Bobroff et ai, Phys. Rev. Lett. 83, 4381 (1999) and 



cond-mat/0010234 . 

P. Mendels et al, Europhys. Lett. 46, 678 (1999). 

A. V. Mahajan et al. Eur. Phys. J. B 13, 457 (2000). 

M.-H. Julien et at, Phys. Rev. Lett. 84, 3422 (2000). 

N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1 694 

(1989); S. Sachdev and M. Vojta, |cond-mat/00092K . 

S. Sachdev, Phil. Trans. Roy. Soc. London, Ser. A 356, 

173 (1998), and references therein. 

N. Nagaosa and P. A. Lee, Phys. Rev. Lett. 79, 3755 

(1997); C. Pepin and P. A. Lee, Phys. Rev. Lett. 81, 

2779 (1998). 

G. KhaliuUin et ai, Phys. Rev. B 56, 11882 (1997); 

R. Kilian et ai, ibid 59, 14432 (1999). 

H. F. Fong et al, Phys. Rev. Lett. 82, 1939 (1999); 

Y. Sidis et al, ibid 84, 5900 (2000). 

S. Sachdev, C. Buragohain and M. Vojta, Science 286, 



2479 (1999); M. Vojta, C. Buragohain and S. Sachdev, 
Phys. Rev. B 61, 15152 (2000). 

D. Withoff and E. Fradkin, Phys. Rev. Lett. 64, 1835 
(1990); L. S. Borkowski and P. J. Hirschfeld, Phys. Rev. 
B 46, 9274 (1992); C. R. Cassanello and E. Fradkin, 
Phys. Rev. B 53, 1507 9 (1996) and 56, 1 1246 (1997); 



K. Ingersent and Q. Si, cond-mat/9810226 

K. Chen and C. Jayaprakash, J. Phys.: Condens. Matter 

7, L491 (1995); R. Bulla, Th. Pruschke, and A. C. Hew- 

son, ibidQ, 10463 (1997). 

K. Ingersent, Phys. Rev. B 54, 11936 (1996); 

C. Gonzalez-Buxton and K. Ingersent, ibid 57, 14254 
(1998). 

N. Read and D. Newns, J. Phys. C 16, 3273 (1983). 

A. C. Hewson, The Kondo Problem to Heavy Fermions, 

Cambridge University Press, Cambridge (1997). 

D. L. Sisson et ai, Phys. Rev. B 61, 3604 (2000). 

We suspect that it is necessary to invoke the coupl ing to 
the dynamic antiferromagnetic spin fluctuations [P0|J 
to robustly explain the observed staggering |l|. 
O. Ujsaghy et ai, Phys. Rev. Lett. 85, 2557 (2000). 
M. I. Salkola, A. V. Balatsky, and J. R. Schrieffer, Phys. 
Rev. B 55, 12648 (1997); H. Tsuchiura et al, Phys. Rev. 
Lett. 84, 3165 (2000). 

E. W. Hudson et ai, preprint. 



